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ABSTRACT 



Context. Classical T Tauri stars (CTTS) are young, late-type objects, that still accrete matter from a circumstellar 

disk. Analytical treatments and numerical simulations predict instabilities of the accretion shock on the stellar surface. 

Aims. We search for variability on timescales below a few minutes in the CTTS TW Hya and AA Tau. 

Methods. TW Hya was observed with SALTICAM on the Southern African Large Telescope (SALT) in narrow-band 

filters around the Balmer jump. The observations were performed in slit mode, which provides a time resolution of 

about 0.1 s. For AA Tau we obtained observations with OPTIMA, a single photon-counting device with even better 

time resolution. 

Results. Small-scale variability typically lasts a few seconds, however, no significant periodicity is detected. We place a 

99% confidence upper limit on the pulsed fraction of the lightcurves. The relative amplitude is below 0.001 for TW Hya 

in the frequency range 0.02-3 Hz in the 340 nm filter and 0.1-3 Hz in the 380 nm filter. The corresponding value for 

AA Tau is an amplitude of 0.005 for 0.02-50 Hz. 

Conclusions. The relevant timescales indicate that shock instabilites should not be seen directly in our optical and UV 

observations, but the predicted oscialltions would induce observable variations in the reddening. We discuss how the 

magnetic field could stabilise the accretion shock. 
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1. Introduction 

; Classical T Tauri stars (CTTS) are young (< 10 Myr), 

low-mass (< 3Mq) pre-main sequence objects. They still 

actively accrete from a circumstellar disk. The disk does 

■ not reach down to the central star, but is truncated by the 

, stellar magnetic field approximately at the co-rotation ra- 

• dius dUchida 1983; Bertout et al. 1988; Kocnigl 1991). The 
[ high-energy photons from the star irradiate, heat and ionise 
, the inner parts of the disk, while the matter is funneled 
' along the magnetic field lines and impacts on the star at 
[ close to free-fall velocity. For the simple case of a dipolar 

• field, the accretion spots are lo cated at high stellar latitude 
I (|Koenig]|ll99ll:IShu et al.lll994{ ). A strong shock develops on 
, or near the stellar surface, which heats up the accreted mat- 
ter to a few 10^ K. Observational support for this model 
comes from different wavelength ba nds. The disk is see n 
in molecular hydrogen emission lines (jHerczeg et al.ll2002t ) , 
and the inner holes can be resolved int erferometrically with 
radio observations (lEisner et al.l 20061). The accre tion fun- 
nel emits strongly in Ha ()Muzerolle et al.l I1998D . This is 
one of the defining characteristics of CTTS. The accretion 
spot on the stellar surface is hot enough to emit UV and 
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X-ray radiation (|Lamzinlll998l:lGiinther et al.ll2007[ ). It also 
heats the un derlying stellar photosphere to about 20 000 K 
(Calvct fc Gullbringj[l998 ). 

This variety of different phenomena leads to time vari- 
ability on different scales. Long changes occur over years, 
possi bly related to a slow ly changing circumstellar extinc- 
tion (jGrankin et al.l 120071 ). For many CTTS the dominant 
source of optical variability is the stellar rotation with 
a typical period of a few days (tRvd gren & Vrba 1983 
iBouvier et allll986t iNguven et alll2009t )'. Depending on its 
inner radius, the rotation timescale of the inner disk and 
also the free-fall time from the disk onto the star is a 
few hours. Thus variability in the accretion rate should 
be limited to this timescale; indeed, variations of about 
0.1 mag are regularly observed in broad-band photo metry 
from CTTS (ISmith et al.lll99fiHRricinski et al.ll2008h . 

Nevertheless, much faster variations can be expected 
if the accretion shock on the stellar surface is unstable. 
This has first been st udied in the context of accretion onto 
white dwarfs (Langer et al.l Il983t IChevalier fc Imamural 
Il982l : Ilmamuralll985t) . Even if the magnetic field is strong 
enough to channel the accretion stream and prevent turbu- 
lence on the boundary between accretion stream and sur- 
rounding photosphere, the shock can still fiuctuate up and 
down. The infalling matter builds up a funnel above the 
stellar surface. Owing to the ram pressure of the infalling 
matter, the gas gets compressed and accumulates. At some 
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point the accumulated matter is dense enough for a run- 
away cooling to set in and the column collapses fast. A new 
shock forms at the stellar surface and works its way upwards 
again buil ding up a new reserv oir of matter. ISacco et all 
()2008() and lOrlando et all ()2009() predict a typical timescale 
of 10 min for the fluctuatio ns, and the same is seen in mod- 
els of lKoldoba et al.l (12008!) , but with much shorter periods 
of 0.02-0.2 s. The difference can be mainly attributed to 
different cooling functions and different input densities n. 
Since the cooling is proportional to n^, diff erent pre-shock 
densit ies heavily influence the cooling time. iKoldoba et al.l 
(1200 8") predict the shock luminosity to fluctuate about 50% 
around the mean, but caution that this value and the exact 
shape of the fluctuation are very model dependent. They 
also find that magnetic fields tend to stabilise the accretion 
shock, so that no oscillations occur. 

The most direct tracer of the accretion shock is the soft 
X-ray emission. This is formed just behind the sho ck front, 
where the accreted materi a l is heated to a 2-3 MK (jLamzinl 
119981: ICxiinther et aHl2007D . iDrake et al.l (|2009h used XMM- 
Newton X-ray data to search for periodicity in TW Hya, one 
of the best-studied CTTS. They did not detect any period 
and set an upper limit on the pulsed fraction of 5% at the 
99% confidence level over the interval 0.001-7 Hz. TW Hya 
is the only target where the X-ray count rate is high enough 
to perform timing analyses on scales below a minute. In this 
paper we use a secondary tracer of the accretion shock. We 
obtained observations in the blue band, using optical tele- 
scopes. In the blue and the UV range, CTTS are domi- 
nated by the continuum emission from the accretion shock 
dLvndcn-Bcll fc Pringle 1974; Grinin 1980: Muzerolle et al.l 
l2000t l because they are late-type objects and their photo- 
spheric emission quickly drops towards shorter wavelengths. 
It is not clear from the current models how fast the opti- 
cal emission reacts to shock instabilities, but observations 
in the visible potentially provide a diagnostic with high 
signal-to-noise ratio, and yield a complementary view on 
sub-minute variability. 

Wc will first introduce the targets observed in Sect. [5] 
and the observations performed in Sect. |31 Details of the 
data reduction process are given in Sect. |31 In Sect. [S]we 
search for periodicities in the data employing Lomb-Scargle 
periodograms. We discuss our results in Sect. [5] and sum- 
marise them in Sect. [71 



2. Targets 

2.1. TWHya 

TW Hya is the closest known CTT S with a distance of 
only 57 pc (|Wichmann et all |1998[ ). Because of the ab- 
sence of a surrounding molecular cloud, TW Hya is a 
key system for the study of CTTS. Br oad Ha profiles 
(FWH M ~ 200 km s~^) were observed bv iMuzerolle et al.l 
(|2000D . TW Hya's stellar paramet ers are M^ ^ 0.7 Mq, 
i?* ~ 1.0 i?0, and its age is 10 Myr ("Webb et al."l999( ). The 
spectral type of T\V Hya is K7 V-Ml V (Webb et alJ [T999t 
iBatalha et al.ll2002[). It is long known that the system is 
seen nearlv pole-on (iKastne r et al. ll997HWilner et al.l2000t 
lAlencar fc Batalhal l2002t ! Weinber ger et al.ll2002D witbjhe 
most precise inclination of 7 ± 1° measured by iQi et al.l 
(|2004 ). Moreover, TW Hya displays variations in line pro- 
files and veiling, which have been interpreted as signa- 
tures of accretion spot rotation (iHerbst fc Koretl 119881 : 



Herbst et all 11994 iMekkadenI 119981: lAlencar fc Batalhal 
20021 : IBatalha et al.ll200i[ 

Ahhough relatively old for a CTTS (« 10 Myr), it is stiU 
actively accreting from a surrounding disk. Modelling of the 
spectral energy distribution from T W Hya requires a gap in 
the inner disk (|Calvet et al.ll2002D and lEisner et al.l (|2006|) 
interferometrica lly resolve its d i sk an d confirm the size of 
the inner hole. iHerczeg et al. (l2002f) observed TW Hya 
in the UV and iHerczeg et al.l (12004) constrain the in- 
ner disk temperature to about 2500 K using H2 fluores- 
cence in HST/STIS and FUSE data. Early photomet- 
ric observations show variability between magnitude 10.9 
and 11.3 in the V-band (Rucinski fc Krautter 1983) . The 
most detailed photometric study so far found variability 
on timescales of hours to days with flicker noise probabil- 
ities, where the amplitude of variations a depends on the 
frequency / as a oc 1/V7 for peri ods between ten day s 
and a few hours (lO^^ - lO""* Hz) (iRucinski et al.l 120081 ). 
TW Hya also was an early targ et for X-ray observations 



with ROSAT (ICosta et al.''20od), Chan dra (IKastner et all 



2002; Raassgaj2009; Brickhouse et al.l |2010D and XMM- 
Newton (IStelzer fc SchmittI l2004ll . In t he X-ray range no 
period was found and lDrake et al.l (l2009( ) placed a 99% con- 
fidence limit that, at most, 5% of the signal is pulsed in the 
frequency range 0.0001-6.81 Hz. 

2.2. AA Tau 

AA Tau is located in the Taurus- Auriga star forming region 
in a distance of about 140 pc. According to its parameters 
it is a typical CTTS with spectral type K7 (1 Cohen fc Kuhil 
Il979l ). a mass of ^ 0.8 Mq, and radius ~ 1.85 Rq. In 
contrast to TW Hya, AA Tau is observed nearly edge-on 
(O'Sulli van et al.l[2005l) and periodically a warp in the disk 
occults the star, causing brightness variations of 1.4 mag 
(iHerbidlTgTTi: iVrba et al.l [19891 [1991 iBouvier et aD [19991: 
iMenard et al.ll2003l ). The period is 8.22 days, which seems 
to be the stellar rotation period ([Bouvier et al.ll2007[ ). The 
occultation is not strictly periodic, but also variable over 
longer timescales. In some instances the fading of the star 
did not occur: this hints at changing accretion geome- 
tries (iBouvier et al.l 12003). The occultation of the UV and 
optical emission coincides with incre ased column densi- 
ties seen in X-ray observations (Schmitt fc Robrada 120071 : 
iGrosso et al.ll2007l) . Apart from one bright, coronal flare 
the X-ray luminosity is not sufficient to repeat the analysis 
done for TW Hya. 



3. Observations 

3.1. TW Hya 

TW Hya was observed with the Southern African Large 
Telescope (SALT) for 45 min in total starting at UTC 23:59 
on 28-03-2009. The spherical SALT primary mirror with a 
total diameter of about 11 m is segmented, fixed in eleva- 
tion, and tilted about 37° from the vertical. This means 
that the telescope can observe only an annulus of the sky, 
which is 12 degrees wide. The CCD camera, SALTICAM, is 
positioned in a prime focus "payload" , which moves across a 
virtual spherical focal surface during an observation, which 
can last up to «2 hours. The full observation of TW Hya 
was scheduled in one crossing. To achieve the highest pos- 
sible time resolution we used the slot mode of SALTICAM, 
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where most of the CCD is covered by a mask and only a 
small strip of 144 rows is read out. With an on-chip bin- 
ning of 6 X 6, an effective sampling period of about 0.14 s 
is reached. SALTICAM consists of two CCD chips with 
two amplifiers each. The CCDs are separated by about 
1.5 mm and each of the amplifiers has a different bias level. 
Observations were taken in three different bands for 15 min 
each. The filters used, in order of observation, are U-band 
and two narrowband UV interference filters with central 
wavelength 340 nm (FWHM 35 nm) and 380 nm (FWHM 
40 nm). TW Hya is severely overexposed in the U-band, 
thus this data was not included in the analysis. At the 
time of the observations SALTICAM operated with soft- 
ware version 4.78. More details o n SALT and the inst ru- 
ment SALTICAM are presented in iBucklev et al.l ()2008D . 

The design of SALT has some implications for photom- 
etry: Only differential photometry is possible since the en- 
trance pupil of the telescopes moves with respect to the 
primary mirror array, thus continuously changing the effec- 
tive collecting area. Also, the light is partially vignetted in 
the slot mode due to the relative position of the slit mask 
and the CCD. Thus the intensity varies as a function of the 
CCD row. 

3.2. AA Tau 

AA Tau was observed with the high speed photo- 
polarimeter OPTIMA Burst (Optical Pulsar TIMing 
Analyzer) on the 1.3 m Cassegrain telescope of the Skinakas 
Observatory (SKO) in Greece on 07-10-2008 from 02:15 to 
03:27 UT. OPTIMA Burst has several possible configura- 
tions. The photometer consists of a single photon counting 
fibre-fed system that contains six fibre apertures which are 
placed in a hexagonal form around the central fiber for 
the target. The apertures have a diameter of 300 microm- 
eters each. This corresponds to 6 arcsec on the sky with 
the SKO telescope. An additional fibre is placed about 1 
arcminute away from the fibre bundle and is used to deter- 
mine the sky background. OPTIMA Burst has an intrinsic 
photon-arrival-time resolution of 4 microseconds and can 
record white light radiation in the wavelength range from 
450 to 900 nm. AA Tau was observed with a blue filter 
(350-550 nm). The fibres are placed in a slanted mirror 
in the prime focus of the telescope. The incoming light is 
also partly reflected by the slanted mirror to a fast-readout 
Apogee AP6 CCD camera with a Kodak chip of KAFIOOOE 
type. It continuously takes exposures with exposure time 
10 seconds while the telescope is tracking. They are used 
to position the target in the central fiber of the hexagonal 
bundle, thus controlling the tracking, and to monitor the 
seeing conditions during the observations. The arrival times 
of the photons in the fibre system are correlated with a GPS 
timing receiver. A more detailed descriptio n of the techni- 
cal pr operties of OPTIMA Burst is given bv lKanbach et ahl 
(|2008l) . 



4. Data reduction 

4.1. TWHya 

SALTICAM data reduction was performed with the 
PySALT packagcQ, which is provided by the SALT 
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Fig. 1. Flux in the filter centred on 340 nm for TW Hya and 
the comparison star. The red/grey line shows a polynomial 
fit. 



Astronomy Operations team and runs in the PYRAF en- 
vironment, a Python extension to the well-known IRAF 
system. Due to a bug in the SALTICAM software up to 
version 4.78 the timing information of the exposures was 
inaccurate. This was corrected using the task slotutcf ix. 
We merged the data read out from all four amplifiers into 
single frames with slotmerge, because the best available 
comparison star is recorded on a different amplifier than 
our target. We extracted the target and comparison star 
using a simple square aperture with side length of 10 pixels 
in each coordinate to cover the entire point-spread function 
(PSF). The background was median filtered and estimated 
from an annular box centred around the target with inner 
and outer dimensions of 16 and 24 pixels, respectively. The 
fluxes of TW Hya and one (340 nm filter) or three com- 
parison stars (380 nm filter) were extracted with the task 
slotphot, which centres the extraction region on the target 
for each exposure. In a few exposures the algorithm cannot 
identify one of the stars and extracts a background region 
only. Those frames can easily be identified because the flux 
drops dramatically in the respective lightcurve for a single 
exposure. Those frames are removed from the lightcurve. 

TW Hya is a factor of 10 and 6.5 brighter than the 
brightest available comparison star in the filters centred 
on 340 nm and 380 nm, respectively. Thus, the errors on 
the relative photometry are dominated by the noise in the 
lightcurve of the comparison star. To detect periodic signals 
we use the higher signal-to-noise lightcurve of TW Hya di- 
rectly and search for periods in the fiux ratio of both stars 
independently. The occurrence of the same period in more 
than one star would be a strong indicator for an instrumen- 
tal effect. 

We show the lightcurves of TW Hya and one compar- 
ison star for both filters in Figs. [1] and [2j The lightcurve 
of TW Hya in Fig. [2] contains two steps after about 300 s 
and 650 s, which are not easily visible in the lightcurve of 
the comparison star. The later step is likely present there, 
too, but hidden in the higher noise level. These steps are 
due to repositioning the target in the slot (an autoguider 
was not yet available). Such a repositioning did not take 
place during the observation in the 340 nm filter. TW Hya 
and the comparison star are affected differently, because 
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400 
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Fig. 2. Flux in the filter centred on 380 nm for TW Hya 
and a brightest of three comparison stars. The red/grey line 
shows a polynomial fit (see text for more details). 
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Fig. 3. Count rate for AA Tan, binned to 1 s (black - left 
axis) and FWHM of a comparison star in the CCD image 
with 10 s time resolution (red/grey - right axis). 



they have different coordinates on the CCD chip, thus the 
steps are still present in the relative lightcurve. This in- 
troduces spurious signals for long periods and we therefore 
limit our analysis to the frequency range 0.02-3 Hz in the 
340 nm band and 0.1-3 Hz in the 380 nm band. We fitted 
the lightcurve of TW Hya with a third order polynomial 
over the entire 340 nm exposure (red/grey line in Fig. [T|) 
and with separate third order polynomials for each step in 
the 380 nm exposure (red/grey line in Fig. [2]) to remove 
trends before further analysis. 

4.2. AA Tau 

We use the data from the central fiber only, which con- 
tains essentially all the signal. The lightcurve shows some 
drop-outs, where no data was read out and a strong in- 
crease towards the end of the night. We selected a clean, 
continuous region of about 40 min before the twilight sets in 
for analysis. This part of the lightcurve is shown in Fig. |31 
binned to 1 s. The second line shows the evolution of the full 



TW Hya - 340 nm 

TW Hya - 380 nm 

comp star - 340 nm 

comp stor - 380 nr 




timescale [s] 

Fig. 4. Autocorrelation coefficient for TW Hya and the 
brightest comparison star in the field in both bands. 

width at half maximum (FWHM) as fitted to a comparison 
star in the accompanying CCD images (timing resolution 
10 s). Times with large FWHM always correspond to low 
flux because light is lost from the instrument. We cannot 
correct for this effect by means of a comparison star, be- 
cause the CCD requires an exposure time of 10 s, which 
is much longer than the sampling rate of the target. The 
lightcurve is mean subtracted before further analysis. 

5. Analysis 

We determine the typical timescale of fluctuations by the 
autocorrelation function. To search for periodicity we flrst 
calculated Lomb-Scargle periodograms of our lightcurves. 
Additionally we performed a wavelet analysis, because pe- 
riodic signals with changing frequency would be washed out 
in a simple periodogram. We quantify the signiflcance of our 
non-detection by analytic considerations and Monte-Carlo 
simulations. 



5.1. Autocorrelation function 

The autocorrelation function is defined as 



P(L) = 



Ejif-^x,x 



i+i 



n 



=N ^2 
= 1 ^3 



(1) 



where Xj are the values of the de-trended and mean- 
subtracted lightcurve and N is the total number of points. 
L is the scale of the correlation. A correlation coefficient 
of 1 indicates a perfect correlation, the absence of a lin- 
ear correlation and -la perfect anti-correlation. We show 
the auto-correlation coefficients for TW Hya in Fig. S] for 
both filters. The autocorrelation of scale is always 1; in 
TW Hya the coefficient differs significantly from the coef- 
ficient calculated for a comparison star for 2-4 s. The au- 
tocorrelation function for AA Tau is different from zero for 
about 10 s, however we flnd similar results using an in- 
trinsically stable white dwarf from a different observation 
as comparison object. Apparently 10 s is the timescale of 
atmospheric variability at the SKO for certain weather con- 
ditions. The lightcurves of the comparison stars for AA Tau 
are not sufficiently densly sampled to test for this. 
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5.2. Lomb-Scargle periodograms 

We calculated normalised Lomb-Scargle periodograms from 
our lightcurves wit h an IDL codqj, which follows the 
method described in lScar glc (1982) and lHorne fc BaliunasI 
(|l986l) . The power P at a given angular frequency uj is cal- 
culated as 



TW Hya - 340 nm 



Pic) 



1 Ej^J COSu(tj-T)^ 

2 y cos^ ^{ij ^'^) 

2 \ sin^ "^(^j —t) 



(2) 



where j runs over all N flux values Xj observed at time tj . 
T is defined by 



tan(2wT) = 



y sin 2ujtj 



I >, cos 2a; ij 



(3) 



In order to normalise the pe riodograms they are scal ed with 
the total observed variance (jHorne fc Baliunaslll986[ ). If the 
measurement errors are Gaussian distributed, the statis- 
tics allow us to calculate a false alarm probability (FAP) 
for each power level. The FAP is the probability that a 
lightcurve of pure noise exhibits at least one frequency with 
the specified power. The value 1-FAP can thus be treated 
as the significance of the detection of any periodic signal. 

For the simple case of a sinusoidal signal with ampli- 
tude A and angular frequency cjq, which is hidden in noise 
with variance a\, , the no rmalised power Pn{u!o) is given by 
(|Horne fc BaliunasllT986[) : 



Pn[ujo) = — 



5.2.1. TW Hya 



9^ 



(4) 



The samphng rate of SALTICAM is 0.14 s, thus the Nyquist 
frequency is just above 3 Hz. In Figs. [5] and [B] the pe- 
riodogram of the lightcurves of TW Hya and the relative 
lightcurve TW Hya and comparison star are shown. The 
large peaks at low frequencies in the second figure are due 
to the steps in the lightcurve, which we could not fully 
correct. No period is detected at the 99% confidence level. 
There are less significant peaks around 0.135 Hz in the peri- 
odogram of the relative lightcurve and 0.14 Hz and 0.155 Hz 
in the absolute lightcurve in the 380 nm filter. The expected 
uncertainty in frequency for the given lengt h of the obser - 
vation is of the order of mHz (our ean.Hland lKovacall981l) . 
so these frequencies differ significantly, where real signals 
should agree; they could possibly be due to aliasing from 
the large low-frequency peaks. In Fig.[7]one periodogram is 
plotte d on logarithmic scales. In contrast to lRucinski et al.1 
(|2008l see their Figs. 3 and 7), who probe periods from 
hours to weeks, we do not observe a flicker noise but find 
the power level to be independent of the frequency. 

5.2.2. AA Tau 

Figure [5] shows a Lomb-Scargle periodogram of the AA Tau 
lightcurve, which was binned to 1 s, and of the correspond- 
ing background data. AA Tau shows significant peaks at 

^ available at 
http : //astro . uni-tuebingen . de/sof tware/idl/aitlib/ 
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Fig. 5. Lomb-Scargle periodogram of the lightcurve of 
TW Hya (black) and the relative lightcurve of target and 
comparison star (red/grey) taken in the 340 nm filter. 
Dashed lines indicate false alarm probabilities. 
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Fig. 6. Lomb-Scargle periodogram of the lightcurve of 
TW Hya (black) and the relative lightcurve of target and 
comparison star (red/grey) taken in the 380 nm filter. 
Dashed lines indicate false alarm probabilities. 
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Fig. 7. Same as Fig. [5] on logarithmic scales 
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Fig. 8. Lomb-Scargle periodogram of the lightcurve of 
AA Tau (black) and the background lightcurve (red/grey) . 
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Fig. 9. Fourier transform of the lightcurve of AA Tau sam- 
pled at 100 Hz. 



very low frequencies, which are also seen in the background. 
There is another, highly significant, peak at 0.038 Hz; how- 
ever, this is a known instrumental signal probably caused 
by feedback between the autoguider and the telescope posi- 
tioning. We conclude that there is no significant periodicity 
between 0.05 and 0.5 Hz. A Lomb-Scargle periodogram of 
the full dataset is computationally too expensive, instead 
we calculated the Fast Fourier transform (FFT) of the full 
lightcurve sampled at 100 Hz (Fig. ^. The largest peak 
seen is again the systematic signal at 0.038 Hz (as already 
discussed) and we can exclude significant periodicity down 
to 50 Hz, which is the Nyquist frequency for a sampling 
rate of 100 Hz. 



5.3. Wavelet analysis 

If periodic oscillations exist in the accretion shocks of CTTS 
their frequency may be unstable. Therefore we performed 
a wavelet analysis of the lightcurves using the IDL func- 
tion wv_applet. Evenly sampled lightcurves are required 
as input, thus missing data points were filled with repeated 
values of the last valid fiux measurement. We employed 



different wavelet shapes, concentrating on Morlets of or- 
der 4-6. For TW Hya no signal with a significance above 
90% is found in the wavelet analysis, while AA Tau shows 
a peak in the power spectrum after about 500 s, which is 
very likely caused by the dip in the lightcurve (Fig. [H]). We 
showed above that this feature is related to changing seeing 
conditions. 



5.4. Upper limits on periodicity 

In TW Hya we do not detect periodic signals on the 99% 
confidence level, that is power level Pjv = 17.4. Knowing 
the variance of the noise aj^ we now invert eqn. 2] to calcu- 
late the minimum amplitude Am which should have been 
detected at that level, if present: 



A-m — CTJV 



/ No 



\APi 



N 



(5) 



For the parameters of both lightcurves from TW Hya Am ~ 
IOctat. The total noise level in the normalised lightcurves 
is about 10^^ and 5 x 10~^ in the 340 nm and 380 nm 
observations of TW Hya, therefore the 99% confidence limit 
on Am is 10~^ and 5 x lO""*, respectively for the presence of 
a single, sinusoidal mode in the frequency range 0.02-3 Hz 
in the 340 nm band and 0.1-3 Hz in the 380 nm band. 

We performed Monte-Carlo simulations to estimate the 
detection efficiency for a combination of many periods. We 
used the time bins from the observations, keeping the irreg- 
ularities introduced by the outlier rejection and simulated 
1000 lightcurves, calculated Lomb-Scargle periodograms for 
each of them, and determined the largest peak in the power 
spectrum. For random measurement uncertainties of 1%, 
which corresponds to the observations of TW Hya in the 
340 nm filter, we distributed multiple sinusoidal waves with 
random frequencies between 0.02 and 3 Hz and added them 
to a constant component. In the case of 40 different frequen- 
cies, each of them with an amplitude of 5% of the constant 
component, more than 90% of the simulations produced 
at least one peak in the power spectrum above the 10% 
FAF, that is higher than in our observations (85% of the 
simulations exceeded 99% FAF). For smaller amplitudes of 
only 2% of the continuum level the peaks in the power spec- 
trum are smaller, still 92% of the simulations with 10 differ- 
ent frequencies produced at least one signal above the 10% 
FAP (91% of the simulations exceeded 99% FAP). Due to 
the larger signal and consequently lower uncertainties these 
values are even better for the data taken with the 380 nm 
fiher. 

Due to the lower signal-to-noise ratio, the data for 
A A Tau is not quite as sensitive. Equation [5] shows, that 
a signal of 0.13 of the measurement error, that is about 
5 X 10"'^ of the total fiux, should have been detected at 
99% confidence level in AA Tau. Since the noise level is 
about three times larger than in TW Hya correspondingly 
larger amplitudes are needed to detect multiple frequencies. 

6. Discussion 

Our analysis of the data presented above yields no signifi- 
cant periodicity in the searched frequency range. The lim- 
its set by this data on the maximum amplitude for varia- 
tions are up to two orders of magnitude more sensitive than 
those iDrake et al.l (|2009r ) derived for the X-ray lightcurve 
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of TW Hya. However, the soft X-rays are formed directly 
in the shock cooling zo ne with a weak contribution from 
the stellar corona (Gii nther et alj l2007f) . the refore much 
stron ger oscillations are expected in X-rays (jSacco et alJ 
|2008[) than in the UV wavelength range we probe. So far 
no time-dependent radiative-transfer calculations of accre- 
tion shocks in CTTS have been carried out. In station- 
ary 1-D models the hot gas at the lower end of the ac- 
cretion column and the heated photosphere below emit 
a sign ificant portion of their e nergy around the B aimer- 
jump (jCalvet fc Gullbrinall998[) . In this wavelength region 
the accretion spot dominates over the ste har photosphere 
by on e order of magnitude in TW Hya (jMuzerolle et al.l 
I2OOOD . The heated photosphere is optically thick and can be 
treated as a black-body w ith temperatures bet ween 5000 K 
and 8000 K (see Fig. 3 in lCalvet fc Gullbrind [l998'). its lu- 
min osity should thus be o f the o rder of 10^^ erg s~^ cm^^. 
The ICalvet fc Gullbrina (|1998f ) model refers to younger 



be seen in absorptio n. In the time-dependent simulations 
of ISacco et al.l (|2008f ) the accretion shock oscillates up and 



CTTS with a radius of 2i?Q, however the spectral type 
of TW Hya matches the model, thus we can apply it to ob- 
tain an order-of-magnitude estimate. The model implies a 
column mass of the heated photosphere down to the region, 
where the continuum is formed, of the order 10 g cm~^, con- 
taining a thermal energy of 10"'^^ erg s^^ cm^^, and suggests 
a timescale for cooling of the heated photosphere of 100 s. 
This calculation gives only a rough estimate, because the 
cooling would change the temperature and thus the opacity 
structure, so the deeper layers become visible as the cool- 
ing proceeds. It suggests that rapid changes of the shock 
conditions will not be visible in optical and UV emission 
from the heated photosphere; however, detailed radiative 
transfer calculations are required to confirm this. 

According to the models the flow time of accreted 
matter through the post-shock zone is a few seconds. 
Incidentally, the flow time through the shock cooling zone 
matches the time we flnd in the autocorrelation function 
for TW Hya. The short term variability in the lightcurves 
shown in figures[T][3]is dominated by statistical fluctuations, 
so the source variability, which is detected in the autocor- 
relation function, is of the order of 1% of the continuum 
level at most. Several effects could be responsible for that: 
If the accretion flow is inhomogeneous and clumpy, this 
would cause the UV emission from the shock cooling zone 
to vary within seconds. The size of the accreted clumps is 
presumably small in comparison to the total area covered 
by accretion spots, which is estimated from 0.3% to 5-10% 
of the stellar surface for T W Hya, depending on the wave- 
length of the observation ([Giinther et al.ll2007l and refer- 
ences therein) . If the clumps switch on or off the emission 
from 1% of the accretion spot, then their characteristic size 
is 10^ — 10'* km. The sound speed in the hot parts of the 
post-shock cooling zone is of the order 100 km s^*, thus 
an area with radius 10^ km could be affected within a few 
seconds; in cooler regions of the post-shock cooling zone, 
where most of the optical emission originates, the tempera- 
ture and thus the sound speed are about an order of magni- 
tude lower. In this region only very small areas can respond 
coherently to disturbances. 

Instead of inhomogeneous accretion the sound waves 
could also be caused by accretion shock instabilities. Three- 
dimensional simulations are needed to study this in more 
detail. But would the optical emission react fast enough 
to changing conditions, however they are caused? Even if 
it does not, the oscillations of the accretion shock could 



down. Thus, the pre-shock column density changes with 
time by ~ 10^° cm~^. This accreted material is not re- 
processed in the hot shock yet and we apply the standard 
conversion between X-ray column density an d optical ab- 
sorpt ion Ay = JVff/2 xlO^* cm~^ = 0.05 mag ("Vuon g et al.l 
I2OO3I ). According to iFitzpatrickl (|l999;) an additional ab- 
sorption of Ay — 0.05 mag would reduce the flux at 3800 A 
by 7%, much more than our upper limit on periodic vari- 
ability. The numbers given are calculated for a line-of-sight 
along the accretion flow. TW Hya is seen pole-on, so this 
condition is valid for accretion at high latitude. The extinc- 
tion changes with the viewing geometry; the argument does 
not apply if the accretion column is seen from the side. 

While general variability is seen, clear periods are ab- 
sent with signiflcant limits down to 0.001 for the maxinmm 
amplitude of a pulsed fraction of the emission. It might 
simply be hard to o bserve oscill a tions in the optical, as 
discussed above, but iDrake et al.l (|2009f ) also do not see a 
periodicity in the X-ray lightcurve of TW Hya. Those au- 
thors present a detailed discussion on the expected period 
of the oscillations, the effect of multiple accretion streams 
and the influence of the magnetic fleld, arguing the sup- 
porting pressure of the magnetic fleld has to be taken into 
account. 

Different studies of time-depen dent accretion shocks 
concentrate on different aspects. iKoldoba et al.l (J2008|) 
present ID magneto-hydrodynamic models of super- 
Alfvenic accretion flows with infall velocities around 
100 km s~^ and different incident angles. In their simula- 
tions strong magnetic fields lead to a stable accretion fiow. 
They develop a stability criterion which depends on the an- 
gle between the field and the stellar surface. However, this 
criterion has not been verified for higher infall velocities. 
If it is still valid for velocities of 525 km s~* and densi- 
ties of 10~*^ g cm~^, as measured from X-ray spectroscop; 



in TW Hya ('Kastner et al.' '2002; 'Stelzcr & Schmitt' "200^^, 
iGiinther et~ l. 2007; Raassen 2009; Brickhousc et al. 2010^, 
then a magnetic field of 100 G is sufficient to keep the shocks 
stable for an accretion flow perpendicular to the stellar sur- 
face. Yet, it is not clear if this is also valid for the stronger 
magnetic flelds usually measured in the accretion s pot re- 
gions of CTTS (jValenti fc Johns-Krulll 12004: Donati et al 
2008), which indicate sub-Alfvenic flows. lOrlando et al. 



(j2010( ) performed 2D simulations of accretion perpendicu- 
lar to the stellar surface including heat conduction in both 
the super-Alfvenic and the sub-Alfvenic regime. They con- 
flrm the importance of the magnetic fleld, and find that 
shocks with high /3 (the ratio between gas pressure and 
magnetic pressure) become s table through magnet ic damp- 
ing for super-Alfvenic flows. lOrlando et al.l (|201(][ ) observe 
quasi-periodic shock oscillations for cases with larger mag- 
netic fields when the flow becomes sub-Alfvenic and chaotic 
behav iour without p e riodic ity for even weaker magnetic 
flelds. lOrlando et al.l (|201(] D also compare their results to 
ID simulations with a similar setup by Sa cco et al.l ()2008[) 
and flnd that 2D models predict a smaller amplitude and 
higher frequency than ID models. For large values of /? 
the magnetic field cannot confine the accretion stream to a 
funnel and matter flows out side- ways, disturbing the sur- 
rounding stell ar atmosphere. 

Al though IKoldoba et all ()2008D and lOrlando et all 
(|2010[ ) predict stability in different regimes, they both 
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agree that there are magnetic field configurations that in- 
hibit sh ock osc illations in the super-Alfvenic regime. Only 
lOrland o et alj ([2010) pr obe sub-Alfven i c flow s, but they 
test less parameters than lKoldoba et all (|2008D . so we can- 
not exclude that similar stability regions exit for sub- 
Alfvenic flows. So, we may not see periodic oscillations of 
the accretion shock, not because of confusion of a larger 
number of accretion streams, but because the shock is in- 
trinsically stable. 



7. Summary 

We presented observations of TW Hya with SALT in two 
filters centred just above and below the Balmer jump with 
a timing resolution of 0.14 s and observations of A A Tau 
with the single-photon counting device OPTIMA. TW Hya 
shows an autocorrelation on timescales of a few seconds, 
however we detect no periodic signal in TW Hya or AA Tau. 
We place a 99% confidence upper limit on the pulsed frac- 
tion of the lightcurves as 0.001 for TW Hya in the frequency 
range 0.02-3 Hz in the 340 nm filter and 0.1-3 Hz in the 
380 nm filter. The corresponding value for AA Tau is an am- 
plitude of 0.005 for 0.02-50 Hz. We simulate lightcurves to 
calculate upper limits for scenarios with multiple indepen- 
dent accretion shocks with different periods. The timescale 
seen in the autocorrelation function is consistent with the 
flow time through the accretion shock and the dimension 
of affected area as estimated from the amplitudes in the 
lightcurve. Our limits on periodicity are so stringent, that 
we discuss the magnetic field needed to stabilise the accre- 
tion shock. 
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